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We provide a fully analytical description of a two-level atom interacting with a broadband light 
field. The problem we present is a typical example of a physical situation which occurs very com- 
monly in practice, but is not straightforwardly solvable. However, we show in this paper that, 
remarkably, only a limited number of elementary calculations are required to treat the problem. 
The presented results are not only highly intuitive, but also allow us to study the nontrivial nonlin- 
ear response of a two-level atom to incident broadband radiation. 



I. INTRODUCTION 



The interaction of a two-level atom and a monochro- 
matic field has always been a very popular topic in quan- 
tum physics. The attraction of this system is owe to the 
relatively modest mathematical tools needed to describe 
the problem, combined with a rich physical behavior, able 
to accurately predict many interesting phenomena such 
as superradiance Q or the Mollow triplet 0. 

The effect of the incident field on the atom's dynamics 
is described by the optical Bloch equations 0,01 ■ This set 
of equations allows one to obtain expressions for the aver- 
age atomic level populations and coherences in the pres- 
ence of an incident field. If the incident light is monochro- 
matic, the optical Bloch equations are easily dealt with. 
For broadband incident light, however, solving the opti- 
cal Bloch equations is in general a far from trivial but 
physically very relevant matter 0, IS 0> IS lU El ■ For 
example, a common experimental situation we are con- 
sidering is an atom or molecule probed by near-resonant 
light, but at the same time irradiated by a broadband 
light source, such as a trapping laser. 

In this paper, we show how to transparently describe 
the interaction of a two-level system and a broadband 
incident wave. An everyday (and practically important) 
example of a broadband field is the field of spontaneous 
photons emitted by an atom. The spectrum correspond- 
ing to spontaneous emission has a Lorentzian frequency 
distribution, and is therefore broadband. Obviously, we 
intuitively expect that if the spectrum of the incident 
field is far detuned compared to the resonance frequency 
of the atom, or very broad compared to the atom's nat- 
ural lincwidth, no significant interaction will take place. 
Likewise, if the incident field has a very narrow frequency 
distribution, we expect that the interaction will resemble 
the one induced by a monochromatic incident wave. In 
this paper, we show that only a few elementary calcula- 
tions are required to quantify the interaction of such a 
broadband field and a two-level atom. We find that not 
just the incident spectrum itself, but in addition the over- 



lap of the incident spectrum and the natural Lorentzian 
emission line of the atom itself determines the interaction 
strength, confirming our intuitive predictions. 

In order to simplify the calculations, we will impose a 
very general restriction on the incident field. More pre- 
cisely, we consider in this paper the class of statistically 
stationary incident fields. The mathematical simplifica- 
tions they allow for, and the fact that many fields encoun- 
tered in practice are statistically stationary, is the reason 
why they are treated in so many textbooks on, e.g., quan- 
tum optics or magnetic resonance. An important prop- 
erty of such fields is that two-time averages {E{t)E{t')) 
of the field E(t) only depend on the time difference t — t' ', 
which implies that field components at different frequen- 
cies are uncorrelated |ll|. 

Finally, we will apply a standard approximation to 
the evolution of the atom itself: we will average out 
the contribution of highly non-resonating terms. For 
a monochromatic incident wave, for example, this ap- 
proximation simply implies that all components oscillat- 
ing at twice the incident wave frequency are neglected. 
The physical meaning of this approximation is that non- 
energy conserving emission and absorption of photons is 
neglected, or in other words, that "virtual photons" are 
not taken into account. For a monochromatic incident 
wave, this approximation is referred to in the literature 
as the Rotating Wave Approximation (RWA). Since we 
deal in this paper with fields which are in general broad- 
band, we will extend the standard approximation found 
in textbooks, and refer to the extension as the General- 
ized Rotating Wave Approximation (GRWA). 

In the following, we will start by formulating the op- 
tical Bloch equations for our configuration. A Fourier 
transform will lead to a better understanding of the be- 
havior of our system in frequency space. Application of 
the GRWA will return steady-state solutions for the op- 
tical Bloch equations. 

To conclude our paper, we will show that the results 
obtained can be straightforwardly used to express the 
nonlinear response of an atom to a broadband field. This 
response, which can be quantified by the atom's dynamic 
polarizability, is of primary importance in studying atom- 
field interactions, since it allows one to gain full insight 
in the internal atomic dynamics. 
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II. THE OPTICAL BLOCH EQUATIONS 



We consider a two-level atom A with lower level a and 
upper level c, separated by an energy difference hw ca . 
The atom interacts with an incident time-dependent real- 
valued field E(i). The dynamics of the system can be 
described in terms of the density matrix a(t). The time 
evolution of a(t) is given by the optical Bloch equations, 
which can be written, in the electric dipole approxima- 
tion, as 



Ona = -I 



& ca &ac J 

ift(t)(cr cc - a aa 
iSl(t)(a cc - a aa 



Face, (la) 



IV CC , (lb) 



lLU ca a ca - -j°ca-, (Id) 



where the notation f(t) = if it) is used. Equations JTJ 
can be found in many elementary books on quantum 
optics [T^ |. The matrix elements a aa and a cc are 
the ensemble-averaged populations of the lower and 
upper atomic level, respectively. They are related by 
o'aa + &cc = 1, expressing conservation of population. 
The off-diagonal elements represent coherences; we will 
elaborate on their physical meaning in section llVl The 
constant decay rate T represents spontaneous emission 
by the system to the surrounding vacuum. The Rabi 
frequency f2(f) = —y t d ac ■ E(t) quantifies the interaction 
strength between the atom and the incident field, 
with d ac the c — > a transition dipole moment. In the 
conventional case of a monochromatic incident wave 
E(t) = Eq cos ujj^t which is often found in the literature, 
the definition £lR a bi = —\d ac ■ Eq is mostly used, 
explicitly removing the oscillatory time dependence of 
the field from the Rabi frequency (obviously, the optical 



Bloch equations then contain extra factors e L ). 
However, we will see that for a more general time 
dependence, as the one we deal with here, it is beneficial 
to use our definition and consider a time-dependent 
Rabi frequency. 

Our aim in this section of the paper is to derive 
steady-state solutions for equations QJ. The statistical 
properties of the field are especially advantageous in 
the frequency domain, since if the two-time average 
(il(t)il(t + t)} of the pha se of only depends on r, one 
can deduce (see, e.g., [llj) that in the frequency domain 



(n[cj]Q[cj']) = J[w]6(u + u'), 



(2) 



where J[uj] is the spectral density function of the incident 
radiation and the Fourier transform of Q(t) is defined as 



1 r +cc 
n[oj] = — / n(t)e- liAjt dt, 

/+oo 
Q,[uo]e lut duo. 
-OQ 



(3a) 
(3b) 



We will see that the appearance of a delta function in 
(J2J) will highly facilitate the calculations further on. In 
what follows, we will Fourier transform the optical Bloch 
equations. We therefore define the population difference 
w(t) and its Fourier transform w[oj] as 



w(t) = a cc (t) - a aa (t) 



w[tu}e- luJt dt, (4) 



and 



f+OO 

&ac{t) = I "acL 1 
J — oo 

/+oo 
a ca [uj]e tuJt duj 
-oo 



-foe 



<7 oc [w]e M *dw, (5a) 



a ca [u]e iwt duj, (5b) 



where the contribution of the non-resonating part of 
the coherences has been omitted. The restriction of the 
integration interval in expressions J5J has the simple 
physical meaning that non-energy conserving terms are 
not taken into account, as mentioned in the introduction 
of our paper. This approximation is a straightforward 
gen eralization of the Rotating Wave Approximation 
[l2| and is therefore referred to as the Generalized 
Rotating Wave Approximation. This generalization 
holds excellently if J[uj] is only appreciably different 
from zero near the atomic resonance. 

If we now split the Fourier transform Q[&] of the Rabi 
frequency Q(t) into a positive- and a negative- frequency 
part 



/■+oo />0 

n(f) = / fl[uj}e iuJt duj + / n[uj]e luJt duj 

JO J -oo 

= n+(t) + n_(t), (6) 



we see that neglecting all highly non-resonant terms in 
results in slightly altered optical Bloch equations: 



(7a) 
(7b) 

(7c) 
(7d) 



CTcc = 


+iQ + (t)cr ca - 


iO_(i)cr oc - Tctcc, 




daa — 


-iVL + (t)<T ca -+ 


iQ-(t)a ac + Ta cc , 




<?ac — 


-iQ, + (t)(a cc 


- cr a a) + iw ca a ac - 


r 

2° ac 
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It is clear that equations are far simpler to deal 
with than the original Bloch equations Q. Fourier 
transforming equation H7a|) and (|7b|) gives 



(r + iuj)w[uj} = +2i / 57[a/]er ca [ti> — ui']dui' 



2i / n[uj'}a ac [uj - lu'}(Lj' -T5(u). 



(8) 

Fourier transforming (j7c|) and (|7d|) . on the other hand, 
gives 



p p + oo 

(iui - iui ca + — )a ac [oj] = -i il[w"]w[w - uj"}duj" , 
z Jo 

(9a) 



(iuj + iuj ca + — )a ca [u>] = +i J H[uj"]w[uj — aj"}du>'' 



(9b) 



If we now substitute the previous expressions in (JHJ), we 
find 



(T + iu)w[u>] +T6(u)) 



/■+00 rO 

2 duj 1 I dcu"n[uj'}n[uj"}w[oj -J - lj"]x 

Jo J — oo 

1 



1 



iuj — iuj' + iuj ca + iuj ~ ilu" — iuj, 



ca T 2 



(10) 



which is a self-consistent equation in the population 
difference. We use the following trial solution for w[oj}: 



w[uj] = w[0}S{oj), 



(11) 



which is appropriate since we are interested in the 
regime for which the populations are time-independent 
(which is also referred to as the "steady-state" regime). 
Substitution of (fTTfl in ((Tn|) and phase-averaging yields 



F 



Tw[0] + T = -2w[0] / dcu'J[uj'} — 



l ( W '- WcQ ) 2 + (|)2' 
(12) 



and therefore 



We can conclude that in steady-state, we find the desired 
solution 



1 - <J„„ = 



X 



2X + 1 

1 

2X + 1 /_ 



(14a) 



-e luJt dcu, 



oo v^ca 



ILUrn — IUJ 



with 



X = 



diuJ[u 



K Q -^) 2 + (^) 2 ' 



(14b) 



(15) 



Expressions l|14f> are the key result of our paper. The 
influence of the spectral properties of the incident field 
enters the dynamics of the density matrix through a sin- 
gle interaction parameter X. The expression (11511 for X 
is surprisingly simple and appealing: it is the overlap in- 
tegral of the spectral density of the incident field, and 
the natural Lorentzian emission line of the two-level sys- 
tem itself. The structure of the interaction parameter 
confirms what we intuitively expect: resonant fields with 
a narrow distribution interact strongly with the atom, 
while the interaction with broad, far off-resonance fields 
is far less pronounced. As examples, we will in the next 
section focus on two specific and interesting values for X. 



III. EXAMPLES 

As a first demonstration of equations l|14(l . we will 
verify that the well-known expressions for an incident 
monochromatic field |l2j | can be retrieved from our 
results. We consider a real-valued monochromatic 
incident field E(t) = E cos(ujLt) with lul > 0. The 
corresponding Rabi frequency is 



(5(u) - uj l ) + 5(lo + Ll> L )^J , 



(16) 



therefore 

fi 2 



(5(oj - uj')6(uj - uj l ) + 8{lo - lo')5{lo + ui L ) 
+ 5(lu + lu')5(lu - uj l ) + 5(lu + lu')S(lu + uj L )j . 



(17) 



Of the 4 terms appearing in (|17f) . only the first remains 
in the (G)RWA, yielding 



W Ld 



i + 2 r 00 dw'j^']—, — 



-6(u). (13) 



J[u] = ^6(u; - lul), 



(18) 



3 of© 



Dynamics of a two-level system in a broadband field 
Tom Savels, Allard P. Mosk and Ad Lagendijk 



which transforms expressions 114() into 



1 



fin 1 



(19a) 



(19b) 

which corresponds exactly to the solutions for incident 
monochromatic fields found in the literature [12| . justi- 
fying our method. 

As a second demonstration of equations l|14l) . we 
consider the nontrivial case of an atom interacting 
with incident spontaneous emission. The incident field 
then has a Lorentzian spectrum with a width T(l + n), 
k > — 1, centered around wr, 3> T: 



JM 



where the factor 



Jo 



§(! + ") 



Jo 



7T ( W - Wi ) 2 + (|(1 + K )) 



J[w]dw = (fi(i) 2 ) 



(20) 



(21) 



to incident light, one has to determine the strength with 
which the atom interacts with the incident light. A key 
property quantifying this strength is the atom's dynamic 
polarizability a. The associated polarization P induced 
by the incident field Eq is then given by P = £n<3 • Eq. 
In the limit of low-intensity incident fields, the response 
is linear and the polarization is then proportional to the 
incident field. However, at higher incident intensities, 
the response is no longer linear since the atom will 
start to saturate. The polarizability will then depend 
on the incident field and the polarization will no longer 
be proportional to the incident field. We will now 
show how the results from the previous sections can be 
straightforwardly used to quantify the rate at which the 
atom saturates due to a incident (broadband) field. In 
other words, we will show how expressions (|14|) allow us 
to calculate the nonlinear dynamic polarizability of a 
two-level atom. 

We take the general case of an incident field E(t) 
consisting of a monochromatic cosine component 
with amplitude Eq and frequency u>l, and a non- 
monochromatic component Ei,(t), with a Lorentzian 
frequency distribution. The fields E(t) thus obtained 
define a very general class, containing the special cases 
of incident monochromatic light, as well as pure incident 
spontaneous emission. The Lorentzian component of the 
incident field has a lincwidth L(1 + k), k > — 1, centered 
around ujl T: 



is proportional to the total incident field energy. We 
find as an explicit solution for X: 



E[u>] = E ^(s(ui + + 5{u - u L fj +E l [u]. 



(23) 



X = 



+ OC 



du}J[u 



1 



-f-oc 



duj- 



i 



i 



(w - u L Y + (L(l + K ))2 ( W _ u; ca )2 + (1 ) 



L^2 



7T I 



Jo 



(2 + «) 



K-c ca ) 2 + (%-) 2 (2 + K )2' 



(22) 



where the extension of the integral in H22fl from to 
—oo is clearly justified by ojl 3> T. Equation 122|) 
clearly shows that the atom-field interaction is weak 
for spectrally broad or far-detuned fields, as mentioned 
earlier. To conclude this paper, we will in the next 
section use the result (|14|) with (|19fl and \\T2\\ obtained 
here to determine the nonlinear response of an atom to 
incident broadband light. 



The corresponding Rabi frequencies are 



—Ml — d ac ■ E , 
-hQ L [u] = d ac ■ E L [u], 



obeying 



(24a) 
(24b) 



{Q L [oj]n L [uj'}) = J l [u]6(lj + u'), 



(25) 



with Ji, [to] the spectral density of the incident Lorentzian 
field 



JlM 



J(k) 



(1 



(26) 



IV. APPLICATION 

Let us now focus on a practical application of the 
previous sections. If one studies the response of an atom 



Equations l(251) - (|2lfl) fully describe the incident field. 
We will now focus on the response of the atom to this 
incident field. Since the dynamic polarizability and the 
density matrix both describe the response of the atom 
to incident light, both quantities should be connected. 
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Indeed, they are related by 



dynamic polarizability of a two-level atom irradiated by 
a monochromatic held is recovered 



e a(u}) ■ E[uj] = d ac a ca [-uj}. uj > 0, (27) 



Our general result l|14|) now allows us to deduce that 
(j23 leads to 



e a(uj) ■ \^-E 5(uj + uj l ) + E l [uj] 



d ac / 
K + l \ 

dgc t 

K + 1 V 



1 Qq5(uj + ujl) 



2X 
d 
2X 



2 uj ca - UJ - ij LOca—U- ij 

1 d ac ■ E q 6(uj - lo l ) d ac -E L [u] 



LUm -UJ — 1- 



UJ rn - UJ - I- 



(28) 



where ® represents the tensor product of two vectors. 
The nonlinear polarizability is thus expressed as 



a(ui) 



d ac ® d ac 1 



e h 2X + 1 
_ 1 



W m - UJ - It; 



Z UJ - UJ ca + in + 



2S 



(29) 



where the static polarizability is given by 



ao 



dac ® ^ac 



(30) 



and where the saturation parameter can be written in a 
surprisingly simple way as 



O 2 
"0 



<S = — +/ dwJ L [w]- r ■ (31) 

4 Jo w-w ca ) 2 + (^) 2 



Expressions l|29|) - (|31|l fully describe the response of a 
two-level atom to an incident field of the general class 

Two limits for the incident field are in particular 
interesting. Firstly, for Jl [u>\ — » 0, the expression for the 



a(u>) = -a ~ 



(32) 



LO — U) ca —lT 



Secondly, for Qq — > 0, we find the saturation due to 
pure spontaneous emission 



a(uj L ) = -S 



2(2X + 1)uj l - UJ ca + «S 



(33) 



with X given by expression (|15|l . 

Finally, we note that expressions (|32|l and 133|) imply 
that for small incident fields, the same expression 



a(uj L ) = -H 



1 



2 UJL - + '• 



(34) 



for the linear dynamic polarizability jl4j is obtained, as 
it should be. 



V. SUMMARY 

In this paper, we have solved the optical Bloch equa- 
tions for a two-level system interacting with a statisti- 
cally stationary broadband field. The resulting steady- 
state density matrix is similar to the result found for 
a monochromatic incident field; the difference between 
both results can be intuitively understood. Finally, we 
have applied the obtained results to calculate the full re- 
sponse of a two-level atom to a broadband field, expressed 
by the atom's dynamic polarizability. 
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